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Using  ray  theory,  a general  analysis  is  presented  for  the  treatment 
of  the  effects  of  sound-source  motion  on  the  total  acoustic  field  at 
a fixed  receiving  point.  Sound  speed  is  depth-dependent  and  the  cw 
source  follows  an  arbitrary  path  *’  arbitrary  velocity.  The 
received  signal  is  interpreted  as  .■>  waveform  of  arbitrary  frequency 
whose  amplitude  and  phase  are  time-dependent.  Application  of  the 
theory  is  made  to  a constant  sound-speed  channel  in  which  the  source 
follows  a short  straight-line  path  with  constant  speed,  and  a 
linearized  acoustic  model  is  developed.  Both  primary  and  cumulative 
acoustic  phase  are  examined  as  functions  of  time  for  various  source- 
trajectory  directions.  An  averaged  Doppler-shift  frequency,  as  well 
as  unaveraged  source  frequency,  is  considered.  Maximum  phase  change 
occurs  for  radial  source  motion  and  minimum  change  occurs  for 
circumferential  motion.  Finally,  acoustic  phase  is  investigated  when 
the  source  follows  a long  straight-line  path.  An  approximate  formula 
is  derived  for  cumulative  phase,  which  is  found  to  be  a hyperbolic 
function  of  time. 


FORM  * 


Ray  Theory 

Underwater  Acoustics 
Moving  Sound  Source 
Acoustic  Phase 


Acoustic  analysis  for  a moving  source  with  application 

to  an  isospeed  ocean 

G.M.  Jacyna  and  M.J.  Jacobson 

Rensselaer  Polytechnic  Institute,  Troy,  New  York  12181 
J.G.  Clark 

Institute  for  Acoustical  Research,  Miami,  Florida  33130 

Using  ray  theory,  a general  analysis  is  presented  for  the 
treatment  cf  the  effects  of  sound-source  motion  on  the  total 
acoustic  field  at  a fixed  receiving  point.  Sound  speed  is 
depth-dependent  and  the  cw  source  follows  an  arbitrary  path 
with  arbitrary  velocity.  The  received  signal  is  interpreted 
as  a waveform  of  arbitrary  frequency  whose  amplitude  and  phase 
are  time-dependent.  Application  of  the  theory  is  made 
to  a constant  sound-speed  channel  in  which  the  source  follows 
a short  straight-line  path  with  constant  speed,  and  a linearized 
acoustic  model  is  developed.  Both  primary  and  cumulative 
acoustic  phase  are  examined  as  functions  of  time  for  various 
source-tra jectory  directions.  An  averaged  Doppler-shif t 
frequency,  as  well  as  unaveraged  source  frequency,  is  considered. 
Maximum  phase  change  occurs  for  radial  source  motion  and  minimum 
change  occurs  for  circumferential  motion.  Finally,  acoustic 
phase  is  investigated  when  the  source  follows  a long  straight- 
line  path.  An  approximate  formula  is  derived  for  cumulative 
phase,  which  is  found  to  be  a hyperbolic  function  of  time. 


INTRODUCTION 


Ray  theory  has  been  used  extensively  in  the  study  of  acoust^  j 

transmissions  in  an  underwater  channel,  where  the  position  of 

both  a point  source  and  point  receiver  are  assumed  fixed.  Aspects 

of  the  moving  source  problem  have  been  treated  in  relation  to  its 

12  3 

effects  on  radar.  ' ' However,  results  are  not  generally 

applicable  in  an  underwater  environment  because  refraction, 
boundary  reflection,  the  coherent  summation  of  multiple  arrivals, 
and  other  effects  are  not  usually  considered.  Only  recently 
have  investigations  been  directed  toward  consideration  of  the 

4 5 

moving  source  problem  in  underwater  acoustics.  ' Computer  results 
for  total- fie  id  phase  and  amplitude  were  presented  in  the  case 
of  a siource  moving  radially  outward  at  fixed  speed  front  a fixed 
receiver,  when  a specific  deep-water  sound-speed  profile  is 
assumed . 

The  intent  of  this  paper  is  twofold.  We  wish  to  present  an 
analytical  technique  whereby  a stationary  sound-source 
problem  can  be  transformed  into  a moving  source  problem, 
ircluding  a careful  explanation  of  Doppler  effects  for  a cw 
sound  source.  In  addition,  to  illustrate  the  technique,  we 
consider  a specific  moving  source  problem  in  an  isospeed  channel. 
The  velocity  of  the  source  is  chosen  constant,  but  its  direction 
with  respect  to  a fixed  receiver  is  not  restricted  to  the  radial 
case.  The  total  acoustic  field  is  analyzed  for  both  short-and 
long-range  runs  of  the  source. 

In  Sec.  I,  a general  ray-theory  treatment  of  the  moving 
source  problem  is  investigated,  in  which  a general  technique 


is  formulated  for  determining  Doppler  effects  for  a cw  source 
moving  with  arbiJ -ary  velocity  over  an  arbitrary  smooth  curve 
in  three-dimensional  space.  A discussion  of  the  effects  of 
source  motion  on  ray  geometry,  boundary  reflection  theories, 
spreading  loss,  and  travel  time  leads 

to  total-field  equations.  In  an  illustration  of  the  general 
approach,  Sec.  II  deals  with  a specific  moving  source  problem 
for  the  case  of  a constant  sound-speed  medium,  bounded  by  a 
horizontal  surface  and  bottom,  where  both  the  sound  source  and 
receiving  point  are  located  on  the  ocean  bottom,  and  where  the 
source  is  assumed  to  follow  a straight-line  path  with  constant 
speed.  General  equations  developed  in  Sec.  I are  used  to  construct 
the  total-field  phase  and  amplitude  as  a function  of  receiver 
time.  In  Sec.  Ill,  certain  restrictions  are  imposed,  under 
which  a linearized  model  of  the  relationship  between  acoustic 
phase  and  receiver  time  is  developed  and  explored  in  detail. 

Phase  variation  for  a source  following  a long  straight-line  path 
with  constant  speed  v is  examined  in  Sec.  IV.  In  particular, 
a relation  between  phase  and  receiver  time  is  developed  and 
used  in  sampling  phar-»  along  the  entire  source  trajectory.  In 
addition,  an  approximate  expression  for  relative  cumulative 
phase  for  all  receiver  times  is  derived.  Section  V summarizes 
the  principal  rerults  of  this  paper. 

I . GENERAL  TREATMENT 

We  begin  by  taking  a receiving  point  (ft  to  lie  at  the 
origin  of  a Cartesian  system  in  which  the  x and  y axes  are 
horizontal  and  the  z axis  is  vertical.  Suppose  that,  at  any 


4 


instant  in  time,  a point  sound  source  5 also  lies  in  the  x-y 
plane  and  has  velocity  with  corresponding  speed  denoted  by  v. 
It  is  well  known  that,  for  "line  of  sight"  transmission,  as 
shown  in  Fig.  1(a),  the  frequency  f detected  at  the  receiver 
differs  from  the  frequency  f emitted  at  the  source.^  If  the 
distance  between  the  source  and  receiver  remains  constant,  the 
wavelength  detected  at  is  >.  = c/f,  where  c is  the  (constant) 

speed  of  sound.  However,  if  the  source  moves  as  shown  in 
Fig.  1(a),  then  at  any  instant  in  time  it  is  effectively  moving 
away  from  CR.  witn  a speed  v cos  , when  05ji><_ir/2  rad.  The  motion 
is  toward  (R  for  r-/2<'j'v  rad.  If  R describes  the  receiver- 

— — /*A 

source  position  vector  at  any  instant  in  time,  where  R is  the 
corresponding  distance,  then  the  resultant  wavelength  detected 
at  the  receiver  a time  R/c  later  is  equal  to 

>R  = (c+v  COSv)  / f , (1) 

and  the  resultant  frequency  is 

fR  “ c/AR  ~ ftl  + (v/c)  cost)-]  . (2) 

In  the  case  of  radial  source  motion,  we  find  that 

fR  = f (1  1 v/c)-1  , (3) 

where  the  plus  sign  corresponds  to  ip  = 0 and  the  minus  sign  to 
“ r . In  essence,  the  term  v cos^  can  be  interpreted  as  the 
scalar  projection  of  the  velocity  v onto  the  negative  ray 

AAA 

direction,  where  in  this  case  the  negative  ray  direction  is 


5 


that  of  R.  Letting  v denote  this  scalar  projection#  we  may 
aav  P 

write 

fR  - f(l  ' vp/c)“l  . (4) 

The  same  technique  can  be  applied  when  a moving  source 
follows  an  arbitrary  path  T with  arbitrary  velocity  in  three- 
dimensional  space  where  the  sound  speed  is  a function  of  depth# 

the  source  and  receiver  are  at  arbitrary  depths,  and  the  source 

/\ 

depth  may  change  with  time.  Let  R(t)  denote  the  position  vector 

AM 

A A ^ 

from  the  receiver  to  the  source  at  a time  t,  and  let  ^v(t)  * djj^dt 

denote  the  corresponding  velocity  vector,  as  shown  in  Fig.  1(b). 

Even  when  sound  speed  is  constant,  the  possibility  of  boundary 

reflections  leads  to  a family  of  rays  which  propagate  from  the 

source  to  the  receiver.  If  we  let  0 (t)  denote  the  angle  at 

n 

the  source  of  the  nth  ray,  measured  from  the  x’  direction,  and 

/V 

let  u (t)  denote  the  unit  vector  associated  with  the  negative 
,w*n  3 

A 

nth-ray  direction  at  time  t as  shown  in  Fig.  1 (b) , then 

A A A 

u (t)  = ccs  0 (t)  i'  + sin  9 (t)  k , (5) 

iwn  n n aaa 


where  i’  and  k are  unit  vectors  in  the  x'  and  z directions, 

/W\  AX'. 

respectively.  Note  that  x*  measures  horizontal  directed 


distance  in  the  vertical  plane  containing  the  ray  and  that 


°n  < ^or  a ra^  directed  upward  (downward)  at  the  source. 

A A 

The  scalar  projection  of  v(t)  onto  u (t)  is 

AAV  AAA»n 


v (t)  = v (t)  • u (t)  = J v ( t ) | cos  tp  (t)  , (6) 

P AV.  MAn  ** 


where  the  angle  p (t)  is  determined  by  v and  u when  their 

MA  AVvP 


j 


6 


initial  points  coincide.  Applying  Eq.  5 to  Eq.  6, 


v (t)  = (v  (t)i+v  (t ) j+v'  (t ) k]  • i cose  (t)i'tsine  (t)k] 
p x m*  y 2 fM  n /<*.  n 


= [R  ( t ) v ( t ) + R ( t ) v ( t ) ] cos0 
ax  y y 


^(t)+Ry(t)] 


-1/2 


+v  (t)sinf  (t)  , 
z n 


(7) 


where  i,  j,  and  k ate  unit  vectors  in  the  x,y,  and  z directions, 

/W\  /V* 

respectively,  and 


i' 

aM 


R (t)i 

A NV\ 


Ry(t) 


R^(t) 


R2  ( 

y 


;>] 


-1/2 


(8) 


is  the  unit  vector  in  the  positive  x'  direction.  A subscript 
on  R and  v refers  to  the  component  of  R and  v in  the  corresponding 
direction.  When  the  source  path  T lies  in  the  x-y  plane, 
v^  = 0 and  Eq.  7 simplifies  to 

O.  A.  .*v  /*\ 

v (t)  = r (t ) • v (t) cosO  (t)  , (9) 

p '•Vi  m«\  n 

where  r(t)  =:  i'  . 


Since  the  source  is  moving,  the  frequency  associated  with 

the  nth  ray  detected  by  the  receiver  at  a time  t is  not  the 

frequency  corresponding  to  the  source  position  at  t.  Rather, 

it  is  that  corresponding  to  the  source  position  at  an  earlier 

time.  Letting  t describe  the  time  associated  with  the 

continuous  emission  of  the  nth  ray  from  the  source,  there 

exists  a relationship  between  source  time  t and  receiver 

r n 

time  t given  by 


t 


n 


t ~ (t  ) , 
n n 


(10) 


- 7 - 


where  Tn  is  the  travel  time  of  the  nth  ray.  Observing  from  Eq.  10 
chat  t is  an  implicit  function  of  t»  we  have 


t * Fn(t)  * 

n n 


(id 


so  that  Tn(tn)  may  be  rewritten  as  T^(t).  Clearly,  the  frequency 
of  the  nth  ray  detected  at  a time  t is  the  frequency  of  the  ray 

-■>  A 

emitted  at  the  time  t * t-T  (t) . Evaluating  Eq.  7 at  t*t  and 

n n n 

using  this  transformation, 

vp(tn)  * vp(t)  " [Rx(t)vx(t)+R  (t)v  (t)]  cos0n(t)  [r^ (t ) +Ry (t )] 

•<•  •» 

> v (t) sinG  (t)  , (12) 

z n 


"1/2 


where  R (t)  * R (F_(t)],  v (t)=v  [F  (t)],  and 

x,y,2  x,y,2  n x,y,2  x,y,z  n 

6 (t)  = 0 [F  (t) ] . The  corresponding  frequency  at  the  receiver, 
n n n 


analogous  to  Eq.  4,  is 
Rn 


f„_(t)  = f [l+{  [Rx(t)vx(t)+R„(t)v..(t)]cos0_(t)  [R?  (t)+R.2,(t))'1/2 


>1 


y y 

-l 


n 


+ v (t)  sine  (t)  }c  (t) 
z n 


y 

(13) 


where  c(t)  represents  the  sound  speed  at  the  source  depth  at 
receiver  time  t.  If  Eq.  9 is  used  in  the  case  where  V lies  in 
the  x-y  plane,  then  Eq.  13  simplifies  to 

-1 


f_  (t) 
Rn 


= f [l+r(t:  -v(t)|  cose  (t)jc_1(t)]  . (14) 


At  the  receiver,  the  contribution  of  the  nth  ray  at  time  t is 


Vn(t)  ‘ VVsin  i“Rn(t»VSnVi 
= An(t)sin  |“Rn(t)[t-Tn(t)J 


+ (t) 

n 


9 


(15) 


where  t = t-T  (t)  is  given  by  Eqs.  xO  and  11,  A (t)=A  [F  (t) ] 
nn  nnn 

is  amplitude,  (t) =S^ [F  (t) j is  phase  change  resulting  from 

surface  and/or  bottom  reflections,  if  any,  and  w (t)=2Trf  (t) . 

th  Rft  Rn 

If  the  n ray  experiences  no  boundary  reflections,  then  Sn(t)H0. 

Substituting  Eq.  13  into  Eq.  15  with  a)=2irf,  we  obtain 

Vn(t)=An  (t)  sin  [u  j l+[Rx(t)vx(t)+Ry(t)vy(t)]c"1(t)cos0n(t) 

r «s.  9 .*  9 -l  -1/2  ^ ^ n | “1  ^ 

x (V  (t)+Ry(t)  J +vz  (t)c"i(t)sin6n(t)|  [t-Tn(t)]  +Sn(t)J  . 

(16) 

We  now  turn  to  a brief  discussion  of  the  phase  term  S 

n 

and  the  amplitude  An  in  Eqs.  15  and  16.  If  Rayleigh  reflection 

7 

theory  is  used,  then,  at  any  one  reflection  of  the  nth  ray 
from  a horizontal  boundary, 

i r 22  i1/2r  2 V*1/2  | 

B exp(ie  ) = (p  /p  ) - (c  /c  ) -cos  0 1-cos ‘0  t 

n n ! m w i_  vr  m wJ  L wj  ) 

i r 2 2 -,1/2  r 2 r1/2  I -1 

■illA,*rAl’CMfl.]  v-c°s  6wJ  j • (17) 

In  Eq.  17,  the  subscript  w refers  to  the  water  while  m refers 

to  the  medium  on  the  opposite  side  of  the  reflecting  boundary 

either  the  ocean  bottom  or  the  atmosphere.  The  quantity  Bn 

is  the  amplitude  loss  at  a reflection  of  the  nth  ray,  e is 

n 

the  shift  in  phase  there,  and  the  symbol  i represents  the 
imaginary  unit.  By  Sno/l’s  law, 

COS0  =[c  /c  (t  )]  cosQ  (t  ) = [c  /c  (t)]  cos 8 (t)  . (18) 

wLw  n nn  w n 


When  Eq.  18  is  substituted  into  Eq.  17,  it  becomes  apparent  that 


amplitude  loss  and  phase  shift  at  a boundary  reflection,  in 

general,  are  time-dependent.  In  addition,  the  number  N of 

bottom  reflections  and  the  number  M of  surface  reflections  of 

the  nth  ray  may  be  time-dependent,  so  that  we  should  write 

N=N  (t  )=N  (t)  and  M=M  (t  ) =M  (t) . If  we  make  the  usual 
n n n n n n 

assumption  of  no  amplitude  loss  and  a Tr-rad  phase  shift  at  a 

surface  reflection,  then  the  contribution  of  all  boundary 

reflections  to  A (t)  is  exp  [n  (t)lnB  (t)l  and  S (t)  is  given  by 

n n n J n 

Sn  (t)  = Nn(t)en(t)+irMn<t)  , (19) 

where  B and  e refer  only  to  a bottom  reflection.  If  the  nth 
n n J 

ray  experiences  no  bottom  reflections  at  time  t,  then  Nn(t)sO; 

if  it  does  not  reflect  from  the  surface  at  time  t,  then  Mn(t)=0. 

Of  course,  other  reflection  models  may  be  considered.  If,  for 

8 

example,  a frequency-dependent  model  is  employed,  then  the 
variation  of  frequency  with  time,  as  described  in  Eq.  13,  must 
be  included. 

When  the  sounu  source  and  ocean  medium  are  stationary, 

when  sound  speed  is  depth-dependent  only,  and  when  the  ocean 

channel  possesses  horizontal  boundaries,  the  contribution  of 

-1/2 

spreading  loss  L to  amplitude  is  L , where  L has  been 
9 l n 

written  as  ' 

L = I R ( sin0  / cos 0 )DR/D9  I.  (20) 

1 x o o 1 

In  Eq.  20,  Ci  is  the  ray  angle  at  the  source,  9X  is  the  ray 
angle  at  the  receiving  point,  and  R is  the  horizontal  distance 
between  source  and  receiving  point.  It  has  been  shown11  that 
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a non-stationary  ocean  requires  corrections  in  the  spreading  loss, 

because  a moving  ocean  causes  an  omhidirectional  source  to 

become  directional  and  because  rays  are  no  longer  normal  to 

wavefronts.  However,  the  corrections  are  negligible  in  an 

isospeed,  uniform-current  medium  when  the  ratio  of  current  speed 

to  sound  speed  is  small.  Since  there  exist  some  analogies 

between  a moving  source  and  a moving  medium  problem,  corrections 

in  Eq.  20  are  expected  here  also.  However,  in  subsequent  sections 

source  speed  will  be  a small  fraction  of  the  speed  of  sound  so 

that  Eq.  20  will  be  assumed  to  be  applicable.  In  using  Eq.  20 

for  the  nth  ray,  R should  be  replaced  by  R (t) =R [F  (t) ] and 

0^  by  0 ( t ) =0  [F  (t) ] , where  F (t)  is  given  by  Eq.  11  and  where 
onnn  n 

R=  (R2+R2 ) 1//2 . Further,  from  Snell's  law,  Isin0  I should  be 
x y 1 x1 

- 2 2~  1/2 

replaced  by  {l-[c  /c(t)]  cos  0 (t' } , where  c is  the  sound 

k n k 

speed  at  the  fixed  receiver.  The  contribution  of  spreading  loss 

to  amplitude  is,  therefore, 


L-i/2(t)  = 
n 


R (t ) l-[cR/c(t)]  cos  0_(t) 


1 1/2,  - 

nH  |coS0„<t> 


-1 


x 3R(t)/36n(t) 


-1/2 


(21) 


Thus  the  amplitude  of  the  nth  ray  at  the  receiver  is 

A (t)  = L 1//2  (t ) exp  [N  (t)  InB  (t)  ] , (22) 

n n n n 

-1/2 

where  Bn  is  given  by  Eqs.  17  and  18,  and  is  given  by  Eq.  21. 

With  the  phase  term  S^  and  amplitude  An  known  from  Eqs.  19 
and  22,  respectively,  the  receiver  arrival  Vn  associated  with 


iswaBsssseMO^iL _^::v  :v ':  “.!  '.J'" 
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the nth  ray  is  known  from  Eq.  15  or  16.  In  order  to  determine 
the  total  acoustic  field  at  the  receiver,  we  elect  to  add  and 
subtract  the  factor  u,t  in  the  argument  of  the  sine  in  Eq.  15. 
The  quantity  uj  is  any  constant,  but  arbitrary,  circular 
frequency.  Then 

V (t)  = A (t)  sin  [u)t+<f>  (t)  ] , (23) 

n n n 

where 

<p  (t)  = oj_  (t)t-wt-w_  ( t ) T (t)+S  (t)  . (24) 

n Rn  Rn  n n 

The  total  field  may  be  expressed  as 


A (t)  sin  [cct+<f>  (t)  ] = £A  (t ) sin  [wt+<J>  (t)  ] , 

n n 


(25) 


where  the  amplitude  A and  phase  <J>  can  be  determined  from  the 
following  system  of  equations: 


A (t)  = 


£An  (t ) sin4> 


-1 


n(t>] 


£a  ( t ) cosi)  (t) 
n n 


sin  4>(t)  = A (t)  £a  (t)sin<t>  (t)  , 

n n 


cos  <t>(t)  = A 1 (t)  Z)An  (t)cosc}>n  (t)  . 


(26) 

(27) 

(28) 


Since  A(t)  and  ^(L)  are  time-dependent,  it  is  apparent  from 
Eq.  25  that  the  total  acoustic  field  at  the  receiver  is  not  a 
sinusoidal  function  of  time.  However,  it  may  be  interpreted 
as  a signal  with  constant  circular  frequency  u>=2nf,  having 
time-dependent  amplitude  A and  phase  <!>.  If  f is  chosen  to  be 
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the  frequency  emitted  by  the  sound  source,  then  u>=w=2Trf  in  Eqs. 
23-28.  A second  procedure,  similar  to  that  of  Ref.  4,  is  to 
average  frequency  over  the  first  N*  significant  ray  arrivals 
and  over  a time  interval  t.  If  we  define 


Aw 


(N’t) 


t ' + T 

t’ 


[ iDRn  (t) -a)]  dt  , 


(29) 


where  t'  is  some  constant  value  of  receiver  time,  then  w should 
be  replaced  by  w+Aid  in  Eqs.  23-28. 

The  above  results  are  too  general  to  predict  the  presence 
of  shadow  and  caustic  boundaries  at  any  instant  of  time.  These 
boundaries  are  important,  however,  since  they  can  be  expected 
to  lead  to  discontinuities  in  Doppler  shift  and  diffraction 
effects.  Once  sound  speed,  receiver  location,-  source  path,  and 
source  velocity  are  selected  for  use  in  a particular  application, 
these  phenomena  can  be  examined. 


II.  THE  ISOSPEED  CHANNEL 


In  some  special  situations,  the  moving  source  problem 
simplifies  greatly.  To  illustrate  the  general  treatment  of 
Sec.  I,  we  consider  here  the  case  of  a constant  sound-speed 
medium  of  depth  H,  bounded  by  a horizontal  surface  and  bottom. 
Eoth  the  sound  source  and  receiving  point  are  located  on  the 
ocean  bottom,  and  the  source  is  assumed  to  follow  a straight- 
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line  path  with  constant  speed  v.  Shadow  and  caustic  boundaries 
cannot  occur  in  an  isospeed  medium.  Suppose,  for  definiteness, 
that  the  source  5 is  located  a distance  R from  the  receiving 
point  <R  at  time  t=0-  Also,  let  the  angle  a in  the  ocean 
bottom  describe  the  direction  of  the  source  trajectory  where 
= (MO).  The  source-receiver  geometry  appears  in  Fig.  2(a), 
while  Fig.  2(b)  shows  the  first  two  rays  from  the  moving  source 
to  the  receiving  point.  Note  that  the  index  n equals  the  number 
of  bottom  reflections  of  the  nth  ray,  and  that  the  source 
location  changes  with  n. 

From  geometric  considerations,  the  distance  between  S 
and  (R  at  t^me  £ is  given  by 


R(t)  = | R (t ) 

Av\ 


ip  a p a I 1/2 

I R + (vt)  +2Rvtcosci| 


(30  ) 


The  scalar  projection  of  onto  the  negative  direction  of  the 
nth  ray  (from  receiver  to  source),  given  by  Eq.  9,  may  be 
written  as 


v 

P 


(t) =r (t) • v (t) cosO 
*v\  ,va  n 


(t)=Bvj[R2 


(t)-R2sin2c.]1/2/R(t)[ 


A 

cos0n (t) , 


(31  ) 


where  0 = sgn  [vt+Rcoscx]  and  where  0^(0  describes  the  angle 

/A 

at  the  source  of  the  nth  ray  at  time  t.  As  discussed  in 

Sec.  I,  the  frequency  of  the  nth  ray  detected  by  the  receiver 

at  time  t is  not  the  frequency  associated  with  the  source 

position  at  t but,  rather,  is  that  corresponding  to  the  source 

position  at  time  t = t-T  (t)  , where  T (t)  = T [F  (t)]  by 

n n n n n 

application  of  Eq.  11.  Under  such  a transformation,  Eq.  31 
with  t=tn  may  be  expressed  as 
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vp  (t)»Bv  j(V(t)-R2sir2ci]1/2/R{t)J  cosO^t) , 


(32) 


where  R(t)=R[F  (t) ] and  6 (t)=0  (F  (t)].  The  corresponding 
n n n n 

frequency  at  the  receiver  is,  analogous  to  Eq.  14, 


1+Bv  | [r2  (t) -R2sin2a]1//2/R  (t)  |c  ^coae^tj'j  , 


(33  ) 


where  c is  now  a constant.  Since  the  travel  time  of  the  nth 
ray  is 


W “ ! R2(tn)  + [2(n+l)H]2  J 


1/2 


-1 


(34). 


(t ) , given  by  Eq.  11,  can  be  determined  explicitly. 
Specifically,  from  Eqs.  10,11,30,  and  34, 

tn=Fn(t)  = | [t+RWc_1cosa]- [T2o+2RMtc_1cosa+M2  (t2-T2Q 

. d2  -2  2 .nl/2  I 2 . -1/2 

+ R c cos  a)J  I (1-M  ) , 


where  M=v/c  is  called  the  Mach  number  and 


T 


no 


2 2 
R /c  + 


[2(n+l)H/c] 


(36  ) 


is  the  travel  time  of  the  nth  ray  if  the  source  were  stationary 
at  (x,y)=(0,R).  The  launch  angle  ®n(t)  may  be  determined  from 

cosen(t;=R(t)  R*(t)  + [2  (n+l)H]  i j , (37) 

~ A 

where  R(t)  can  be  found  from  Eqs,  30  and  35  with  t=t  . Thus, 

n 

from  Eq.  33 , the  corresponding  circular  frequency  at  the 
receiver  is 


15 


u0r,  (t)  = ^Ht2  +M2F2  (t)4-2RM(cosa)F  (t)/c  l1/2j  [T2  4-M2F2(t) 

Rn  no  n n j ) L no  n 

1/2  2 22 

+ 2RM(cosa)Fn(t)/c]  ' 4-BM  [R  (cosa)  /c+M  Fn  (t ) 

4 2RM{cosa)Fn(t)/c]  1/2  1 ] , (38) 

where  ( t ) is  given  in  Eq.  35. 

We  now  make  the  usual  assumption  of  no  amplitude  loss 
and  a ir-rad  phase  shift  at  a surface  reflection.  If  Rayleigh 
reflection  theory  is  used  at  each  reflection  of  the  nth  ray 
from  the  horizontal  bottom,  then  the  use  of  Eqs.  18  and  37 
enables  us  to  write  Eq.  17  as 


B (t) exp [ie  (t) 
n n 

['vV-! 

(cm/cw)2-[R(t)/2(n4-l)H]2 

2 ' 

1/2] 

[ ( ft  / p ) 4-  ' 

2 2 

x [ 1 — ( c / c ) J | 

J 

: (cm /c  ) -[R(t)/2(n4-l)H]': 

w m 1 

? 

m w 1 

-1 

m w 

X [1-  (C  /c  ) ] 

w m 

1 J 

• 

(39) 

Although  this  result  could  be  used  in  Eqs.  19  and  22,  we 
choose  to  make  instead  the  simpler  assumptions  of  a 2-dB  loss  and 
a rr-rad  phase  shift  at  each  bottom  reflection.  Then,  the 
contribution  of  all  boundary  reflections  to  the  ampli iude 
An  of  each  ray  arrival  is  10~  * , and  the  phase  Sn,  given 

by  Eq.  19,  becomes 

S (t)  = nTt4-  (n+1 ) tt=  tt  (mod  2 tt ) . (40) 

n 
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From  Eqs.  30,  35,  and  3?,  R and  e^are  known  and  the  contribution 
of  spreading  loss  to  amplitude  is,  from  Eq.  21, 

L~1/2(t)  = c“1|T2o+M2F2(t)+2RM(coscx)Fn(t)/c  j , (41) 

where  F(t)  is  given  by  Eq.  35.  Thus  the  amplitude  of  the 
nth  arrival,  given  by  Eq.  22,  becomes 

A (t)  = io~°‘ln  L'1/2 (t) . (42) 

n n 

with  S , A , tu_  , and  T known,  the  nth  arrival  V , given 
n n Rn  n n * 

by  Eqs.  2 3 and  2 4,  is 

*“1/2 

V (t)=10~°,ln  c-1  |t2  +M2F2  (t)+2RM(cosa)F  (t)/c! 

n I no  n n I 

x sin  [ut+<J>  (t)  ] , (43  ) 

n 

where 

* (t)=ojpT2  +M2F2  (t)+2RM(cosa)F  (t ) /c  |1/2  j [t2  +M2F2  (t) 
n 1 1 no  n n ) | no  n 

1/2  22  22 

+2RM (cosa)F  (t)/c]  + 3M[R  cos  a/c+M  F (t ) +2RM (cosa) 

n n 

*F  (t)/c]1/2[  1 1 1-  It2  +M2F2(t)+2RM(cosa)F  (t) /c! 1/2 1-wt+r . 

n iJlino  n n)  J 

(44) 

The  total  field  is  now  known  from  Eq.  25 , and  its  amplitude  and 
phase  from  Eqs.  26-28. 

III.  LINEARIZED  MODEL 

Under  certain  assumptions,  it  is  possible  to  expand  Eqs.  43 
and  44  in  powers  of  the  Mach  number  M.  In  particular,  an 
expansion  that  keeps  only  linear  terms  in  M greatly  simplifies 
these  equations  while  preserving  the  effect  of  source  motion. 
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When  the  rauio  of  source  speed  to  sound  speed  is  small  (M<<1)  , 
such  a linearization  is  valid  when  restrictions  are  placed  on 
time  duration, and#  subsequently,  source  movement. 

If  one  performs  a linearization  in  M of  Eq.  44 , then 

d>  (t)z  2Mtt  w (cosa)T  -2Mn  u (cosa)  t-wT  + (w-w)  t+tr  , (45) 

n n no  n no 

where 

i 2 

n » R/T  c * l+[2(n+l)H/R]  I . (46 ) 

n no  ( 1 

•** 

A similar  linearization  of  A (t) , given  in  Ec  42,  leads  to 

n 

A (t)  = Kf°*xn  R"1  n [l.-Mn  (cosa)  (t/T  “1)  j . (47) 

n n n no 

A further  simplification  of  Eqs.  45  and  47  is  possible  if  H<<R 

and  n is  not  too  large.  In  such  a case,  n = 1-6  and  T z (R/c)  (1+6) 

n no 

2 

wh^re  6=2  [ (n+1)  H/R]  <<1.  If  nn  is  replaced  by  unity  and  rj?no 
by  R/c  only  in  those  terms  containing  M in  Eqs.  45  and  47,  we 
obtain 

<+>  (t)~2ML)Pc  ^cosa-2Moj  (cosa)  t-coT  + (w-w)  t+ir  , (48) 

n no 

and 

A (t)=  10'°*ln  r'1  n U-M(cosa)  (ct/R-x)  ] . (49) 

n n 

Of  course,  the  remaining  nn  and  TnQ  terms  may  be  approximated 
also,  if  desired.  The  corresponding  total  acoustic  field, 
given  in  Eq.  25,  may  now  be  expressed  as  the  imaginary  part  of 
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A (t ) exp  i(£t+<t  :t)  ] ~ (l-M(coaa)  (ct/R-1)  lexp 
*cosa-2Mu)  (cosa)  t+  (oj-ui)  t ] 10  ^*^nR  ^nn 


i (Jt+2MwRc~* 

expi(-«T  +tt). 
no 


(50) 


Although  the  summation  in  Eq.  25  contains  an  infinity  of 
terms,  its  rapid  convergence  permits  its  termination.  In 
Eq.  50  we  have  kept  N’  terms  where,  for  example,  N’  miaht 
be  chosen  so  that  the  amplitudes  of  subsequent  arrivals 
are  less  than  0.01  that  of  the  n * 0 ray.  The  termination 
justifies  the  earlier  assumption  that  6 is  small  when  H<<R. 
The  summation  in  Eq.  50  represents  the  approximate  total  field 
for  a stationary  source  at  (x,y)=(0,R).  Since  this  location 
corresponds  to  t=o»  we  shall  write  the  corresponding  amplitude 
and  phase  as  AQ(0)  and  <J>Q(0),  respectively.  Thus, 


NL-1 


A_(0)exp  i (0) 
o o 


10  0,ln  r expi  (-wT  +tt)  . 

n no 


(51) 


The  reason  for  our  selective  approximations  in  Eqs.  48  and  49 
is  now  clear:  we  are  able  to  express  the  complex  form  of  the 

total  field  as  the  product  of  a source-motion  factor  and  the  total 
field  for  a stationary  source,  where  each  term  in  the  latter  is 
exact.  If  a term  expiwt  is  cancelled  on  each  side  of  Eq.  50 
and  if  Eq.  51  is  employed,  a comparison  of  the  moduli  and 
arguments  in  the  resulting  equation  leads  to  the  following 
linearized  approximations  for  the  total-field  phase  and 
amplitude : 
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$(t):  4>0 ( 0 ) •*  2Mu(coaa) (t-R/c)+(w-  w)t  , (52) 

and 

A(t):  Aq(0)  [l-M(coso)  (ct/R-1)]  • (53) 

Since  u is  usually  large,  Mu  Is  not  small  even  for  a slowly 
moving  source.  Thus  Eq.  52  suggests  that  the  effect  of  source 
motion  on  phase  variation  is  significant,  in  qeneral.  However, 
the  effect  on  amplitude  variation  is  negligible  if  Mc/R<<1. 

If  u is  chosen  to  be  the  cixcular  frequency  u at  the 
source,  then  Eq.  52  simplifies  to 

♦ (t.)  « $.(t)=  <t  (0)-2Mu  {coset)  (t-R/c)  . (54) 

1 o 

A second  procedure,  outlined  in  Sec.  I,  is  to  average  frequency 
over  the  first  N'  significant  ray  arrivals  and  over  a time 

interval  t.  When  we  carry  out  an  average  analogous  to  Eq.  29, 

we  obtain 

Aw~  - Mw(cosa)  , (55) 

where  the  right-hand  side  is  just  the  D >ppltr  shift  associated 

with  the  n=0  ray.  If  we  take  w=w+Au,  with  Au  given  by  Eq.  55, 

then  Eq.  52  gives 

<t(t)  = ^(t)  ~ CO)  - 2Mu  (cosa)  (t/2-R/c)  . (56) 

We  observe  from  Eqs.  54  and  56  that  both  and  are 
approximately  linear  in  time,  where  we  have  assumed  that  |t| 
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is  small.  Furthermore,  the  total-field  phase  is  the  same  for 
a = aQ  and  a **  - a , so  that  it  is  necessary  to  examine  phase 
only  for  0<a<_n  rad.  Differentiation  of  Eq.  54  with  respect 
to  time  leads  to  the  observation 


(t) 


= - *i(t: 

ot*ct 


a=7r-a 


o 


(57) 


Thus,  the  slopes  of  the  4>,  versus  time  curves  for  a=a  and 

l o 

a=r-aQ  are  equal  in  magnitude  but  opposite  in  sign.  In  other 
words,  source  trajectories  with  a component  approaching 
toward  or  receding  from  the  receiving  point  at  the  same  rate 
give  rise  to  the  same  magnitude  of  phase-rate.  These  comments 
apply  to  as  well.  If  Eq.  56  is  differentiated  with  respect 
to  receiver  time  and  compared  with  the  derivative  of  Eq.  54, 
we  see  that 


'J^(t)  ~ 24-2  (t)  . (58) 

Therefore,  the  rate  of  change  of  phase  with  averaged  frequency, 
v,,  is  only  half  that  of  the  phase  with  unaveraged  frequency, 

V «.  • 

1 

The  relative  phases  and  ^""^o  aPPear  versus  time  in 

Figs.  3 and  4 for  selected  values  of  the  source-trajectory 
angle  u in  deg.  For  these  figures  we  selected  the  parameter 
values  R=10  km,  c = 1500  msec  \ v = 10  kn,  and  f - 350  hz. 
Although  the  water  depth  H does  not  appear  in  Eqs.  52-58,  it 
will  be  recalled  that  our  results  assume  that  H/R  is  small. 


With  the  above  parameter  values  and  H * 100m,  Eqs.  52,  54,  and 

5e  were  found  to  be  in  excellent  agreement  with  the  exact  phase 

given  by  Eqs.  26-28  and  Eqs.  43-44,  when  the  time  interval  was 

restricted  to  approximately  -5.0<t55.Q  min.  However,  Figs.  3 

and  4 show  a time  interval  of  half  this  amount.  The  lower 

horizontal  axes  indicate  the  directed  distance  r in  Km  traveled 

by  the  source  from  its  position  at  time  t * 0.  Figure  3 

illustrates  the  linearity  of  phase  over  a relatively  short 

time  interval  and  the  fact  that  the  slopes  of  the  phase  curves 

for  a=a  ° and  a=(180-a  )°  are  equal  in  magnitude  but  opposite 
o o 

in  sign.  Note  in  Fig.  3 that,  over  the  5.0  min  interval 
shown,  the  phase  variation  is  as  much  as  730  cycles,  where 
the  maximum  variation  occurs  when  the  sound  source  follows 
a radial  path  (a- 0° , 180° ) . When  a»90°,  the  source  trajectory 
is  tangent  to  the  circumferential  direction  at  t*=0  and  the 
approximate  phase  variation  is  zero  as  expected,  since  the 
source-receiver  range  is  virtually  unchanged.  Figure  4 is 
similiar  to  Fig.  3,  but  shows  the  relative  phase  w*ven 

average  Doppler  shift  is  included.  In  agreement  with  Eq. 
curves  in  Fig.  4 have  only  half  the  slope  of  corresponding 
curves  in  Fig.  3. 

Primary  phase  is  phase  modulo  one  cycle,  and  will  be 

indicated  in  this  paper  by  the  subscript  p.  Both  the  primary 

relative  phases  ] and  [$>  -4>  ] appear  in  Fig.  5 over 

lop  2 o p 

the  time  interval  -5 . 0<t<5 . Osec  when  the  source  moves  radially 
outward  from  the  fixed  receiver  (a=0°).  As  in  Figs.  3 and  4, 
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it  is  observed  that  phase  variation  is  decreased  by  one-half 
upon  inclusion  of  an  averaged  Doppler-shif t frequency.  The 
two  curves  are  qualitatively  similiar  to  corresponding  results 
in  Ref.  4,  even  though  phase  variation  in  that  reference  was 
calculated  for  a deep-water  sound-speed  profile  and  for 
different  parameter  values. 

IV.  THE  LONG  TRAJECTORY 

Under  the  assumptions  of  the  previous  section,  we  now 

consider  the  source  to  follow  a long  straight-line  path  with 

constant  speed  v.  We  let  d denote  the  distance  of  closest 

approach  of  source  to  receiver  and  let  t=0  describe  the  time 

at  which  S is  the  distance  d from  <R  . Suppose  now  that 

the  source  is  a distance  R(t  ) from  the  receiver  where  a(t  ) 

o o 

is  the  corresponding  angle  subtended  by  the  source-path  and 
receiver-source  directions,  as  shown  in  Fig.  6.  In  Secs.  II 
and  III,  R represented  the  initial  distance  between  source  and 
receiving  point  at  t=0,  and  <i>o(0)  described  the  phase  of  a 
stationary  source  at  a distance  R from  the  receiver  at 
time  t= 0 in  Eq.  54.  In  that  equation,  w was  chosen  to  be  the 
circular  frequency  at  the  source,  and  we  shall  restrict  our- 
selves to  this  case  here.  The  analog  of  Eq.  54  , for  4>  expanded 

about  t rather  than  zero,  is 
o 

<J>(t)  1 (tQ)  - 2Mcu  [cosa  (tQ)  ] [ (t-to)-R(tQ)/c]  , (59) 

where  |t-t  |<<1  and  4>  (t  ) is  the  phase  of  a stationary  source 
'o’  00 

at  time  t=t  . Now,  let  t be  the  time  required  for  6 to  travel 
0 0 
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from  the  point  P of  closest  approach  to  to  the  generic 

position  shown  in  Fig.  6.  If  the  distance  traveled  is  r and 

if  the  speed  of  8 is  v,  then  r=vt  . Both  R(t  ) and  a(t  ) 

o o o 

are  found  from  the  following  expressions: 

R(to)  = (d2+v2t2]1/2  = c[  (d/c)2+M2t2]1/2  , (60) 

and 

cosa (tQ) = vto[d2+v2t2]"1/2  = MtQ(  (d/c)2+M2t2]"1/2  . 

( 61) 

Using  Eqs . 60  and  61/  Eq.  59  may  now  be  expressed  as 

*(t)  : $o(to)+2M2toto-2M2wtot(d/c)2+M2t2]"1/2(t-to)  (62) 

A 

for  |t-t  I < 5 min  and  H/R<<1.  If  we  let  $ (t,  t ) =4>  (t)  -$  (t  ) 
o'—  o o o 

denote  the  phase  relative  to  that  a stationary  source  at  the 
position  corresponding  to  t=tQ,  then  Eq.  62  can  be  rewritten 
as 

“1/2 

i(t,to)  : -2M20)to  j C(d/c)2+M:t2]  (t-t0)“lj  . (63) 

A 

In  Fig.  6,  primary  phase  [4>(t,t  )]  is  sampled  for  seven 

^ IT 

20  sec  intervals  over  a period  cf  2 hours  or  a source  travel  of 
approximately  37  km,  for  d = 10  km.  When  tQc0,  it  is  seen 
that  the  effect  of  a slowly  moving  source  on  phase  variation 
is  negligible.  Such  a result  is  consistent  for  a source  path 
tangent  to  the  circumferential  direction,  as  discussed  in  the 
previous  section  and  shown  by  the  a=90°  curve  in  Fig.  3.  As 
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| t | increases  from  zero,  a steady  increase  in  the  magnitude 
of  phase  variation  is  observed.  For  large  values  of  |t  |, 
the  source  is  effectively  moving  along  a radial  path  from  the 
fixed  receiving  point,  where  maximum  phase  variation  was 
predicted  in  Fig.  3 by  the  a=0°  and  a=180°  curves. 

A continuous  representation  of  approximate  relative 
cumulative  phase,  rather  than  primary  phase,  can  be  formulated 
since  it  is  reasonable  to  assume  that  all  time  derivatives  of 
the  phase  exist  in  some  neighborhood  of  each  value  of  tQ. 

Then,  to  within  the  order  of  the  approximations  made  in  the 
previous  section,  Eq.  62  represents  the  first  two  terms  of 
the  Taylor  series  for  <Mt)  about  the  point  t=tQ.  Since  the 

I 

Taylor  coefficient  of  (t-tQ)  is  4>  (tQ)  , we  have  from  Eq.  62 
that 

$'(t  ) = - 2M2ut  [ (d/c)2+M2t2]"1/2  . (64) 

o o o 

Integration  of  this  result  with  respect  to  tQ  gives 

*(tQ)  = -2w[ (d/c) 2+M2t2] 1/2+C,  (65) 

where  C is  a constant  of  integration.  Because  we  are 

interested  only  in  relative  cumulative  phase,  we  select 

C = 2wd/c  in  order  to  make  $(0)  = 0.  In  addition,  we  write 

t in  place  of  t to  obtain 
c o 

*(t)  =-2w[  (d/c)2+M2t2]1/2  + 2wd/c  . (66) 


For  large  values  of  |t|,  Eq.  66  becomes 
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^(t)  ~ -2Mw  1 1 j + 2tud/c  . (57) 

Thus  the  magnitude  of  the  rate  of  phase  variation  is 
approximately  constant  for  large  |t|,  as  observed  in  Eq.  62 
or  in  Fig.  3 as  a approaches  0°  or  180°.  Further,  Eq.  67 
may  be  construed  as  an  upper  bound  on  the  cumulative  phase. 

The  approximate  cumulative  phase,  given  by  Eq.  66 , is 
plotted  as  the  solid  curve  in  Fig.  7 over  an  eight-hour  time 
period  or  a corresponding  source-travel  distance  of  148  km. 

The  phase  is  observed  to  be  hyperbolic  and  its  asymptotes, 
given  by  Eq.  67/  are  graphed  as  dashed  lines.  The  figure 
reiterates  previous  conclusions.  For  t ~ 0,  the  source  is 
near  its  point  of  closest  approach  to  the  receiver,  and  phase 
variation  differs  little  from  zero.  As  |t|  increases,  a 
rapid  increase  in  variation  is  observed,  which  approaches  a 
constant  when  the  source  trajectory  is  effectively  radial 
relative  to  the  receiving  point.  The  cumulative  phase,  Eq.  66/ 
is  a global  result  that  was  obtained  from  a local  result,  Eq.  62. 
Since  the  latter  is  only  approximate,  Eqs.  66  and  67  need  not  be 
accurate  for  very  large  values  of  |t|. 

V.  SUMMARY 

This  paper  presents  an  analysis  for  the  treatment 
of  the  effects  of  sound-source  motion  on  the  total  acoustic 
field  at  a fixed  receiving  point  in  the  ocean.  The  source 
emits  a cw  signal,  and  may  follow  an  arbitrary  path  in  three- 
dimensional  space  with  arbitrary  velocity.  Sound  speed  is 
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taken  to  be  depth-dependent.  Using  ray  theory,  it  is  explained 
how  source  motion  affects  ray  geometry,  frequency,  travel  time 
and  spreading  loss,  and  boundary  loss  and  phase  shift.  It  is 
shown  how  the  movement  of  the  sound  source  leads  to  a received 
signal  that  may  be  interpreted  as  a waveform  of  arbitrary 
frequency  whose  amplitude  and  phase  are  t ime-depondent . 

Application  of  the  general  theory  is  made  to  a constant 
sound-speed  (channel  with  horizontal  boundaries  in  which  the 
sound  source  and  receiving  point  are  located  on  the  ocean 
bottom.  The  source  is  assumed  to  follov;  a straight-line 
path  with  constant  speed.  A simplified  linearized  acoustic 
model  is  developed  when  source  speed  is  much  less  than  sound 
speed  and  when  the  ocean  depth  is  much  smaller  than  the  source- 
receiver  range.  The  results  are  valid  over  a sufficiently  short 
time  interval,  and  the  effect  of  source  motion  on  acoustic  phase 
variation  is  shown  to  be  highly  significant.  Both  primary  and 
cumulative  phase  are  examined  as  functions  of  time  for  various 
directions  of  the  source  trajectory.  The  same  magnitude  of 
phase-rate  is  obtained  for  source  trajectories  with  a component 
approaching  toward  or  receding  from  the  receiver  at  the  same 
rate.  Maximum  phase  change  occurs  for  radial  source  motion, 
while  minimum  change  occurs  for  circumferential  motion.  For 
selected  parameter  values,  acoustic  phase  is  shown  to  change  by 
Hundreds  of  cycles  over  a time  interval  of  only  a few  minutes. 
When  the  arbitrary  frequency,  discussed  above,  is  taken  to  be 
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the  average  Doppler-shif t frequency  associated  with  the  first 
ray  arrival,  the  phase  variation  at  the  receiver  is  only  one- 
half  the  variation  when  the  frequency  is  chosen  to  be  that  of 
the  signal  emitted  at  the  source. 

The  above  phase  results,  which  are  local  in  time,  are 
extended  to  the  global  case  in  which  the  sound  source  follows 
a long  straight-line  path  at  constant  speed.  When  the  source 
is  near  its  point  of  closest  approach  to  the  receiver,  phase 
variation  is  small.  However,  a steady  increase  in  variation 
is  observed  as  the  source  departs  from  this  point.  An 
approximate  formula  for  cumulative  phase  is  developed  with 
the  result  that  this  phase  is  a hyperbolic  function  of  time. 

The  asymptotes  of  the  hyperbola  provide  a bound  for  the 
approximate  cumulative  phase  and  an  estimate  to  it  when  the 
sound-source  location  is  far  from  its  point  of  closest  approach 
to  the  receiver. 
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figure  LEGENDS 


FIG.  1. 
FIG.  2. 
FIG.  3. 

FIG.  4. 

FIG.  5. 

FIG.  6. 


Doppler  effect  and  geometry  for  a moving  source  and 
fixed  receiver. 

Geometry  for  linear,  constant-speed  source  motion  and 
two  typical  ray  paths. 

Relative  phase  versus  receiver  time  and  source  travel 
for  selected  source  angles:  w=t d,  R=10km,  c=15oomsec~^ , 
v=10kn,  and  f=350hz. 

Relative  phase  versus  receiver  time  and  source  travel 
for  selected  source  angles:  w=u>+Aw,  R=10km,  c=1500msec 
v=10kn,  and  f=350hz. 

Primary  relative  phase  versus  receiver  time  and  source 
travel  for  the  source  angle  a=0°  and  w»w  and  u=io+A<jj, 
respectively.  Other  parameters  as  in  Fig.  3. 

Primary  relative  phase  versus  receiver  time  and  source 
travel  sampled  over  a 2 hr  interval  for  w=u>  and  d=10  km 
Other  parameters  as  in  Fig.  3. 

Cumulative  relative  phase  versus  time  and  source  travel 
over  an  8 hr  interval  for  w=w  and  d=10  km.  Other 
parameters  as  in  Fig.  3. 
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